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ABSTRACT 

Motivated by the problem of the formation of active regions from a deep-seated solar magnetic field, we con- 
sider the nonlinear three-dimensional evolution of magnetic buoyancy instabilities resulting from a smoothly 
stratified horizontal magnetic field. By exploring the case for which the instability is continuously driven we 
have identified a new mechanism for the formation of concentrations of magnetic flux. 
Subject headings: Sun: interior — Sun: magnetic fields — instabilities — MHD 



1. INTRODUCTION 

Solar active regions are the surface manifestations of a 
deep-seated, predominantly toroidal magnetic field. Although 
it is generally agreed that the solar magnetic field is main- 
tained by some sort of hydromagnetic dynamo, the mecha- 
nism by which this is effected is far from understood. Con- 
sequently, the strength, structure, and location of the interior 
magnetic field are similarly unknown. However, most recent 
solar dynamo models, despite their significant differences in 
other respects, postulate that the tachocline, the thin layer of 
strong velocity shear located at the base of the convection 
zone, plays an important role in the generation of toroidal field 
by the shearing of a weaker poloidal component (see, e.g., the 
review by Tobias & Weissll2007|) . Given this premise, it is im- 
portant to address the nature of the initial escape of the mag- 
netic field from the tachocline, its subsequent ascent through 
the convection zone, and its eventual emergence at the photo- 
sphere. 

Owing to the vast range of scales across the convection 
zone, it is impossible to model realistically all of these stages 
in one calculation. Here we concentrate solely on the insta- 
bility of a layer of magnetic field, with the aim of clarifying 
the physics of the formation of coherent magnetic structures 
from a much larger scale field. 

A vertically stratified, horizontal magnetic field can be un- 
stable to magnetic buoyancy instability provided that the field 
decreases sufficiently rapidly with height. It is important to 
note the non-trivial distinction between instability to two- 
dimensional modes in which the field lines remain straight 
(interchange modes) and instability to fully three-dimensional 
modes. The former are essentially destabilized by a decrease 
with height of B/p (where B is the magnetic field strength 
and p the density), the latter simply by a decrease with height 
of B. The phy sics underlying this difference is elucidated in 
iHughes & Cattaneo (1987). 

The linear theory of magnetic buoyancy instabilities has 
been intensively studied oyer a number of ye ars (e.g. 
lNewcombllT96TI: lParkedlT966l: iGilmanll 19701: lAches on 1978). 
The nonlinear development of the instability, particularly with 
respect to a dee p-seated solar field, has, in evitably, received 
less attention. ICattaneo & HughesI (Il988|) investigated the 
nonlinear evolution of interchange modes resulting from a 
slab of uniform horizontal magnetic field embedded in a con- 
vectively stable and otherwise field-free atmosphere. Here the 
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instability is driven by a density jump at the upper interface 
of the magnetic slab — an extreme form of magnetic buoy- 
ancy instability. This interfacial instability generates a strong 
shear flow that leads, by a secondary Kelvin-Helmholtz in- 
stability, to the formation of strong vortices. The subsequent 
nonlinear evolution is then governed by pairwise vortex in- 
teractions, which can even act so as to drag down pockets of 
strong magnetic field. The three-dimensional instability and 
nonlinear evolution of the same b asic stat e was examined by 
Matthews et al. (1995) and Wissin k et aD ([2000). The initial 
evolution again takes the form of interchange modes, lead- 
ing to the formation of vortex tubes. However, neighbouring 
vortex tu bes of oppos ite sign are unstable to a longitudinal in- 
stability (ICrowlll970l ): this in turn causes the magnetic field 
to adopt an arched structure. By contrast. Fan ( 2001) investi- 
gated the nonlinear evolution of a smoothly varying magnetic 
field, with a profile chosen such that the initial state is unstable 
to three-dimensional modes, but stable to interchange modes. 
She shows the formation of arched magnetic structures, which 
maintain a reasonable degree of coherence as they rise. A 
recent review of the implications of magnetic buoyancy in- 
stabilities for the solar tachocline is given by Hughes (2007|) . 
Related simulations of magnetic buoy ancy in sub- surface re- 
gions have also been performed (e.g. lArchontis et al.l [20041: 
llsobe et al.l[2005[) . 

Inspection of active regions on the solar surface suggests 
the emergence through the photosphere of a buckled, toroidal 
magnetic field. However, very little is known about the struc- 
ture, strength and orientation of the sub-surface field. For 
example, which properties of the observed surface field re- 
sult from the initial instability and which result from later in- 
teractions in the convection zone? Our aim in this Letter is 
to seek further understanding of the instability by examining 
the nonlinear evolution from a very simple equilibrium state. 
We consider an atmosphere that is in both magnetohydro static 
and thermal equilibrium and that has a linear magnetic field 
profile. We concentrate on modes that are intrinsically three- 
dimensional; i.e. the field gradients considered are such that 
the basic state is unstable to three-dimensional modes but sta- 
ble to interchanges. 

Furthermore, we impose boundary conditions such that the 
source of the instability is maintained for all times, thus al- 
lowing us to consider the long-term evolution. This is an im- 
portant difference from the earlier works cited above, all of 
which consider "run-down" experiments, in which the poten- 
tial energy stored in the initial field configuration is rapidly 
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Fig. 1 . — Model with a = ^ and eigenfunction perturbation, (a) Evolution 
of the kinetic energy density. (5) Horizontal average of the magnetic field 
Bx vs. depth, initially (dotted line), at f ~ 50 (dashed line), and in the final 
steady state (solid line). 

converted into kinetic energy and which is then slowly dissi- 
pated. 

We consider four cases, characterised by two different field 
strengths and two different choices of initial conditions. We 
have identified a completely new mechanism for the forma- 
tion of flux concentrations; this is of potential significance in 
the solar context in which the formation of localised regions 
of strong field from a weaker larger scale field is a crucial 
component of the flux emergence process. 

2. MATHEMATICAL FORMULATION AND 
PARAMETERS 

We consider a Cartesian layer of perfect gas governed by 
the equations of three-dimensional compressible MHD. The 
dimensionless representation of these equations is obtained by 
scaling lengths with the layer depth, time with the isothermal 
sound crossing time at the top surface (z = 0), and tempera- 
ture, density, and magnetic field with their respective values 
at the top surface. The physical properties of the fluid are en- 
tirely determined by four dimensionless numbers: the Prandtl 
number Pr, the ratio of viscous to thermal diffusivity; (, the 
ratio of magnetic to thermal diffusivity; 7, the ratio of specific 
heats; and Q, the nondimensional thermal conductivity. 

The basic state that we consider is an exact equilibrium, 
dependent only on depth z, with a unidirectional horizontal 
magnetic field Bx(z) =l-\-z/Hb. The density and temperature 
are determined from solution of the magnetohydro static equa- 
tions. In the absence of a magnetic field, the atmosphere takes 
the form of a plane polytrope, with polytropic index m and 
constant temperature gradient 6. The strength of the magnetic 
field is changed by varying a = 2/P(xB^, where the plasma-/? 
is the ratio of the thermal to the magnetic pressure. A more 
detailed descriptio n of the governing equations and parame- 
ters is contained in lTobias et al.l (|1998|) . 



We consider perturbations to this basic state, subject to the 
following boundary conditions: all variables are taken to be 
periodic in both horizontal directions; the top and bottom 
boundaries are stress-free, impermeable, and isothermal. For 
the magnetic field we choose to maintain a fixed gradient of 
Bx (dBx/dz = l/Hb) at z = 0, 1, together with dBy/dz = and 
Bz = 0. Our aim in this Letter is to examine the nonlinear de- 
velopment of modes that are intrinsically three-dimensional. 
By an extensive linear analysis we have identified the regions 
in parameter space where three-dimensional modes (kx ^ 0, 
ky T^'O) are unstable but interchange modes (kx = 0) are stable. 
In the next section we discuss the nonlinear evolution result- 
ing from two choices of parameters and two different types of 
perturbation to the basic state. 

3. NONLINEAR EVOLUTION OF 
THREE-DIMENSIONAL MODES 

The nonlinear evolution of three-dimensional magnetic 
buoyancy instabilities is studied by solving numerically the 
equations of compressible MHD with a modified version of 
the hybrid finite-dif ference-pseudospectral parallel code used 
by Brumme ll et al.l (1996) (hydrodynamic) and Tobias e t al.l 
([199 8|) (MHD). We consider the evolution from two basic 
states, distinguished by the value of the field strength; the 
other parameters are held fixed at Pr = ( = 2 x 10~^, Ck = 
2.5 X 10-^ = 2, m= L6, Ht = 1, and 7 = 5/3. Further- 
more, we consider two different perturbations to each basic 
state: one consisting of the most unstable eigenfunction, the 
other of random noise. In each case, the horizontal dimen- 
sions of the Cartesian computational box are determined by 
the wavenumbers kx and ky of the most unstable mode such 
that one wavelength fits in the x-direction, along the initial 
magnetic field, and four wavelengths fit transverse to the ini- 
tial field, in the j-direction. 

3.L Evolution to a Steady State 

Here we consider a basic state with a= 1/4 perturbed by 
the most unstable eigenfunction, namely, that with growth 
rate a = 0.16 and with wavenumbers kx = 0.^7 and ky = ^3^ 
(hence, the dimensions of the computational domain are Lx = 
7.24 and Ly = 3). The evolution is essentially linear up to 
time t ^ 21. For 21 < ^ < 37 the kinetic energy continues 
to grow (see Fig.[Tt). During this stage, broad upflows carry 
strong magnetic fields from the bottom of the layer, while nar- 
row plumes pull down weaker field, leading to the formation 
of arched structures, as shown in Figured Figure [T]b shows 
the mean result of this nonlinear evolution, with a significant 
reduction of the gradient of the horizontally averaged mag- 
netic field Bx. Initially constant throughout the layer (dotted 
line), the vertical gradient of Bx becomes nearly zero in a large 
fraction of the domain at ^ ~ 50 (dashed line). The nonlin- 
ear reorganization has therefore acted to remove the driving 
mechanism of the instability. There is an obvious analogy 
with thermal convection, in which saturation is achieved by a 
nonlinear reorganization leading to an adiabatic core with thin 
thermal boundary layers. 

The arched configuration shown in Figure [2l is however 
transient; following this initial nonlinear reorganization the 
system relaxes to a state with weak flows and a magnetic field 
profile that is close to critical. As can be seen from Figure[T]b, 
the field gradient in the final state is reduced everywhere in the 
layer compared with its initial value. This subsequent evo- 
lution is permitted by the choice of boundary conditions, in 
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Fig. 2. — Isosurface of magnetic energy (5^/2 = 1) at f ~ 45 for the same 
model as in Fig.[T] 

which the gradient of Bx (but not the field itself) is fixed at the 
boundaries. 

3.2. Evolution to Time -dependent Concentrated Flux States 

First we consider the evolution from the same basic state 
as in ED (i.e. with a= 1/4), but with the perturbation tak- 
ing the form of small- amplitude random noise. As shown in 
Figure [3^, the linear evolution (up io t ^ ?>6) and the initial 
nonlinear phase, as measured by the temporal evolution of 
gross properties such as kinetic energy, are similar to those 
portrayed in Figure [T^ . The first nonlinear reorganization 
again leads to the formation of arched structures, as shown in 
Figure m although, as a result of the initial conditions, these 
are less regular than those of Figure [2l Crucially, these irreg- 
ularities play a significant role in modifying the subsequent 
behaviour. From 100 < ^ < 450 the kinetic energy increases 
monotonically until the onset of a secondary oscillatory in- 
stability. From t ^ 1200 the system evolves on two disparate 
timescales: a short cycle with period ^ 7 modulated on larger 
timescales 250 < ^mod < 300. 

The underlying mechanism of the short period oscillations 
can be understood by inspection of Figure [3]b, which shows 
the temporal evolution for fixed values of x and z of the mag- 
netic energy, the transverse horizontal velocity v, and the ver- 
tical velocity w. The left panel shows the periodic formation 
of a concentration of magnetic energy, which drifts slowly 
with velocity 0.03. The central and right panels indicate that 
the concentrations of magnetic energy result from convergent 
downflows associated with rolls in the jz-plane. The magnetic 
flux becomes concentrated between the two counter-rotating 
rolls; it thus becomes buoyant and rises rapidly. This drives 
a countercell, which diverges at the flux concentration and 
thereby destroys it. After some reorganization, the initial cel- 
lular flow is reestablished but displaced relative to its initial 
position. Reformed downflows again lead to a concentration 
of flux and the entire process is repeated. 

Next we examine the evolution from an equilibrium with 
a stronger field (a = 1). For this case the evolution is qual- 
itatively similar regardless of the nature of the initial pertur- 
bation; here we shall describe the case in which the pertur- 
bation takes the form of random noise. The most unstable 
mode has wavenumbers kx = 0.S6 and ky = 9J5, and growth 
rate a = 0.45. Figure [5^ shows that the temporal evolution 
is broadly similar to that of the weaker field case portrayed 
in Figure [3^. The secondary instability is evident, leading 
to a modulated periodic state. The increase in initial field 
strength leads, as expected, to a more vigorous instability 
(higher growth rate and increased saturation level) and to a 
shorter period oscillation of the nonlinear state. As in the ex- 
ample above, flux concentration occurs via converging down- 
flows. The peak field in the flux concentrations is significantly 
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Fig. 3. — Model with a = ^ and random perturbation, (a) Evolution of the 
kinetic energy density, (b) Space-time plots, with y horizontal and t vertical, 
of magnetic energy (1.17 < 5^/2 < 1.78), |v| < 0.05, and |w| <0.06; 1863 < 
r < 1913; 0.7. 




Fig. 4. — Isosurface of magnetic energy (5^/2 = 1) at r 74 for the same 
model as in Fig. [3] 

Stronger than the average initial field. In this case, as can be 
seen in Figure [6l the flux concentrations are never completely 
dispersed, although their strength varies as they propagate. 
This (modulated) travelling wave behaviour for the field and 
flows can be clearly seen in Figure [5]b. 

4. DISCUSSION 

In this Letter we have identified a new, inherently nonlinear, 
mechanism for the formation of coherent magnetic structures 
from a layer of weaker magnetic field. The initial develop- 
ment of the magnetic buoyancy instability drives flows that act 
so as to form isolated concentrations of magnetic flux. This 
nonlinear state persists because the instability is continually 
driven from the boundaries and saturation of the instability 
occurs via the net redistribution of flux so as to remove aver- 
age magnetic field gradients from the interior of the compu- 
tational domain. An intriguing feature of the flux concentra- 
tions is that, once established, they travel — taking the form 
of a modulated wave. 

The processes in the deep interior of the Sun that lead to 
the formation of active regions are still poorly understood. 
Even the most optimistic estimates from dynamo theory sug- 
gest that an upper bound for the large-scale dynamo-generated 
field is of the order of 10^ G (see, e.g. Ossendrijver 2003). The 
energy of such a field is small, however, in comparison with 
that of the strong downflows in the convection zone and, such 
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Fig. 5.— Same as Fig.[3] but with o; = 1. In panel b 0.65 < < 2.27, 
|v| < 0.51, \w\ < 0.56; 153 < f < 200; and z - 0.5. ^ ^ 



a field m ight be expected to be pinned down by the convection 
(see, e.g. lTobias et al.lfT998|) . Clearly, therefore, some process 
must be acting so as to create strong localised field structures 
from the dynamo field that are capable of traversing the con- 
vection zone in order to form active regions. Our calculations, 
although idealized, suggest a new mechanism that may play a 
key role in this process. 

We conclude by noting that the discovery of the mechanism 
arose from considering a system where the instability is con- 
tinually driven, in contrast to nearly all previous studies in 
which the initial magnetic buoyancy instability was investi- 
gated by means of run-down calculations. Both types of ex- 
periment may be of importance for understanding the dynam- 
ics of the solar interior. The magnetic buoyancy instability is a 
fast process that might be captured by run-down calculations 
(which usually lead to the formation of arched structures). On 
the other hand, the replenishment of toroidal field by the dy- 
namo should lead to a continual source of instability of field 
at the base of the solar convection zone. 





Fig. 6. — Slices in the yz-plane of magnetic energy at times t ^ 152.9, 
154.7, 156.5. The peak magnetic energy ^ 3. 
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